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I. INTRODUCTION 

Generalized parton distributions (GPDs) Q, |^ |^ ^ 
II II im H Ho"! appear in the context of the QCD fac- 
torization in various hard exclusive phenomena includ- 
ing deeply virtual Compton scattering and hard exclusive 
meson production. Among several general constraints on 
GPDs an important role is played by the polynomiality 
of the Mellin moments |3 and by the positivity bounds 

[illillillillillilEIIlliliSilil. In this 

paper we suggest a representation for GPDs which auto- 
matically satisfies both positivity and polynomiality con- 
straints. The analysis is restricted to the case of spin-0 
hadrons (e.g. pions) but various types of partons will be 
covered. 

We use the following definition of GPDs: 



2^ 



exp(iAx)(P2|0W(A,n)|Pi). (1) 



Here \Pk) is the hadron state with momentum P^,. The 
light-like vector n, 

= 0, 

is normalized by the condition 

n{Pi+P2)^2. (2) 

We use the standard notation of Ji |Q| for the parameters 
A, t and ^ 



A = P2-Pi, e = -^("A), t 



(3) 



The definitions of light-ray operators O^^^ (A, n) for var- 
ious types of partons are listed in the Table HJ We have 
included the scalar field into this table since the pos- 
itivity bounds are more general than their applications 
in QCD. The last column of this table contains the num- 
ber TV of factors n'^ appearing in the light-ray operator 
0''^\\,n). This number N plays an important role in 
the formulation of the positivity bounds and of the poly- 
nomiality conditions and we include N in the notation 
dU of GPD H(^\x,S^,t). 

The structure of the paper is as follows. Section ^ 
contains a brief description of the polynomiality and pos- 
itivity constraints on GPDs. In Section ITTll we introduce 



Parton 




N 


scalar 







quark 


iVi(-^)(n-7)V'(^) 


1 


gluon 


n>^G% (-^) npG^'''P (^) 
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TABLE I: Light-ray operators 0'''^\\,n) for various types of 
partons and the corresponding parameter TV. 



a modified version of the double distribution representa- 
tion for GPDs which slightly differs from the standard 
one but is more relevant for our aims. In Section llVI an 
ansatz for the double distribution is suggested. In the re- 
maining part of the paper we show that this ansatz leads 
to GPDs which obey both polynomiality and positivity 
constraints. In Section we compute the GPDs corre- 
sponding to our ansatz for double distributions. Section 
IVII contains the check of the positivity of the forward 
parton distributions. The positivity bounds on GPDs 
are verified in Section IVIII Appendix ^ contains the 
derivation of the general solution of the positivity bounds 
(without the polynomiality constraints). 



II. POLYNOMIALITY AND POSITIVITY 

Whatever limited our knowledge about GPDs is, there 
are two basic constraints: polynomiality and positivity. 
The polynomiality means that Mellin moments in x of 
GPD ijW(a;,C,t) 



± 

j dxX™ilW(x,e,i) =P,„+Ar(^,t) 



(4) 



must be polynomials in ^ of degree m -\- N . 

Various i neq ualities for GPDs suggest ed in the Refs. 
m m m IM EI EI EI EI E HfliU can be consid- 
ered as particular cases of the general positivity bound 
on GPDs derived in Ref. [23|. This general positivity 
bound has a relatively simple formulation in the impact 
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parameter representation IT^ I20ll2ll[2^ . The impact 
parameter appears via the Fourier transformation of 
the A-'- dependence of GPDs. If the transverse plane 
is orthogonal to vectors n and Pi + P2, then the trans- 
verse component of the momentum transfer A (PJ is 
connected with the variable i = A^ by the relation 



.-L|2 



t 



1-e 



(5) 



Here M is the mass of the hadron. We define the GPD 
in the impact parameter representation as follows: 

{x,tb^)= J |^exp[z(A^6^)] 



lA^ 



1-e 



(6) 



Here notation is used in order to avoid confusion 

with the nucleon GPD H and to keep the compatibil- 
ity with the notation of Ref. [22, where the following 
inequality was derived: 



1 



1-^ 



1 



1+e 



> 0. 



(7) 



This inequality was obtained in Ref. |22| for the case 
iV = 1 and the generalization to arbitrary N is straight- 
forward. 

Inequality ^ should hold for any function p{z). 
Therefore we actually deal with an infinite set of pos- 
itivity bounds on the GPD. The general inequalit y lt7 | 
covers various inequalities suggested for GPDs [ll|, [13, 
ElQIIllIllIlllllllliSilas particular cases with 
some special choice of functions p{z). 

It is well known that the double distribution represen- 
tation ill, 1^ JJ with the D-term j2^ 



dad[35{x - - I3)Fd (a, t) 



|a| + |/3|<l 



+ em-\x\)D 



-,t] sign(0 



(8) 



guarantees the polynomiality property Q. Another in- 
teresting parametrization for GPDs supporting the poly- 
nomiality was suggested in Ref. (25| . 

The positivity bound on GPDs |7J| is equivalent to the 
following representation for GPDs in the impact parame- 
ter representation (see Appendix^) in the region x > |^|: 



JV+l 



1+r 



1-r 



(9) 

with arbitrary real functions Qn- Instead of the discrete 
summation over n one can use the integration over con- 
tinuous parameters. 

Although both polynomiality and positivity are basic 
properties that must hold in any reasonable model of 
GPDs usually the model building community meets a 
dilemma: one can use the double distribution represen- 
tation ^ but it does not guarantee that the infinite set 
of inequalities (0) will be satisfied |2^ . Alternatively 
one can build the models based on the representation 
or on the so called overlap representation 15] which also 
automatically guarantees positivity bounds but then one 
meets problems with the polynomiality. In this paper a 
representation for GPDs is suggested which guarantees 
both positivity and polynomiality. 



III. MODIFIED DOUBLE DISTRIBUTION 
REPRESENTATION 

For the construction of GPDs H^^\x,£,,t) obeying 
both polynomiality and positivity constraints we use the 
double distribution representation which differs from the 
standard representation © by the extra factor of (1- 



H^''\x,U) = il-^) 



N 



X / dad(3S{x-^a- f3)F'jy{a,f3,t). (10) 

q| + |/3|<1 

Here N depends on the type of the parton distribution 
according to Table n 

Representation H1U|) obviously satisfies the polynomi- 
ality condition Indeed, 



da; a;" J dadp5{x - S,a - (3)F'jj{a, (3,t) 

1 |a| + |/3|<l 



= Pn{tt) (11) 

where Pn{S,, t) is a polynomial of degree n. Therefore 

/ 

-1 -1 |q| + |/3|<1 

X 5{x -ia- P)F'o{a, f3, t) = SN+m{L t) (12) 

is a polynomial in ^ of degree N + m in agreement with 
Eq. ©. 
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For our aims it is convenient to use parameters 



ai 



a2 = i(l -/? + «). (13) 



instead of a, (3. Actually it is ai and a2 (oli = 1 — x — y, 
a2 = y in terms of variables x,y used in refs. 0, Qj) 
that appear as a parameters in the perturbative dia- 
grammatic justification of the double distribution rep- 
resentation. The modified double distribution expressed 
in terms of parameters ai, a2 will be denoted as follows: 



FD{ai,a2,t) = F'jj{a,(3,t) . 



(14) 



After these changes the modified double distribution rep- 
resentation H1U|I takes the following form: 



the positivity constraint. We use the following ansatz for 
the modified double distributions H15|) : 

oo 

Fd (ai , q;2 , t) = j dXX j dv 



\a1a2 



t\ L^{Xai,Xa2). (18) 



Our double distribution is parametrized by an infinite 
set of functions L^{wi,W2) defined for wi,'W2 > and 
depending on parameter v. We assume that for any 
v function L^{wi,W2) corresponds to a positive definite 
quadratic form in wi,W2, i.e. for any function (f>{w) 



1 l-Qi 

^ I I da-2 



ii-(^)(2;,C,i) = 2(l-a;)'" / da^ 





X Fd {ai,a2^t)5[x - ^(02 - "i) - (1 - ai 



"2)] 



(15) 

Here we use the triangle integration region in the ai, a2 
plane which corresponds to the constraint /3 > in terms 
of variables a, (3. Hence our GPD vanishes in the "anti- 
quark" region (for brevity we use the word "quark" for 
any type of partons): 



il(^) (x,C,i) = if x<-\£,\. 



(16) 



Therefore we must take care about the positivity con- 
straints only in the "quark" region x > |^|. Once this 
pure quark GPD is constructed we can use the transfor- 
mation X — > —X to build GPDs with appropriate proper- 
ties in both quark and antiquark regions. 

In the case iV > one can add the Z3-term to the 
modified double distribution representation (|15|l : 

i/W (x,e,t)^-^°ifW ix,^,t) 



+ x 



N- 



1 



sign(0 . (17) 



In principle, using the trick of Ref. |23|, one can include 
the i'-term into the double distribution. But if one is in- 
terested in a parametrization of GPDs obeying the posi- 
tivity and polynomiality constraints, then the Z3-term is 
useful: it is localized in the region < |^| and therefore 
it does not appear in the positivity condition |(7J). The 
polynomiality is obvious for the _D-tcrm. Thus by adding 
an arbitrary D-term we violate neither polynomiality nor 
positivity. 



IV. ANSATZ FOR DOUBLE DISTRIBUTIONS 

Now the problem is to find double distributions 
Fo {ai,a2,t) which lead to GPDs TJ'^^^ {x,^,t) obeying 



dwi J dW2Liy{'Wi,W2)4>{wi)(j3*{w2) >0 ■ (19) 



This is equivalent to the existence of the following inte- 
gral representation for L^{wi,W2) 

L„{wi,'W2) =^ J dpF„{wi,p)F* {'W2,p) (20) 

or to its discrete series analog. Since we are interested 
in real and ^-even GPDs, we must use real functions L^, 
and F^. 

The lower limit of the integral over v on the right- 
hand side (RHS) of Eq. |[TH|) determines the asymptotics 
of Fjy {ai,a2,t) at large \t\. If one integrates over posi- 
tive I/, then Fd ^ \t\^^- Functions Li, appearing in Eq. 
(|18|l have the ly dependent dimension, which is slightly 
awkward but simplifies the equations. 

Since X,ai,a2 > and t < the following factor ap- 
pearing in Eq. I|18|) is always positive: 



1 



Xaia2 



-t>0. 



(21) 



Below it will be shown that for any set of positive defi- 
nite functions L^{wi, W2) under the assumption that the 
integrals on the RHS of l|18|) are convergent, the resulting 
double distribution Fd (ai,Q;2,t) (|18|l leads to the GPD 
H^^'^ [x, t) l|15|) which satisfies the positivity bound iQ. 
This check of positivity will be done in Section rVIII but 
first we prefer to derive some useful relations. 



V. EXPRESSION FOR GPDS 

Let us derive the expressions for GPDs H'^'^^ {x,^,t) 
corresponding to the double distribution ifTI^ . First 
we insert ansatz H18|l for the double distribu- 
tion FD{ai,a2,t) into representation (|15|) for GPD 



1 — ai 



7jW(x,^,t) = 2(1-2;)^ / dai 



dao 



4 



x6 [x - ^(a2 - ai) - (1 - ai - 012) 



VI. FORWARD DISTRIBUTION 



X I d\\ I dv 



1 



\a1a2 



-u-l 



-t] L^{\ai,\a2). (22) 



We can rewrite this as follows: 



00 CXD 00 



=2(1-2;)^ d\\ dai da 







x9{l — ai — a2)S [x — £,{a2 — ai) — (1 — ai — a2)] 



X diy 



1 



t] Ly(Aai,Aa2 



(23) 



^ AaiQ;2 

Let us introduce new integration variables 

Wk = Xak (24) 
instead of ak and integrate over A using the delta function 

{x,^,t)=2{l-x)^-' 



00 00 



X / dwi / dw20 { X — ^ 



, W2 - Wi 
' Wl + W2 







In the forward limit ^ — > 0, i ^ we obtain from Eq. 
(El 



/(x)-ijW (x,0,0) = 2(l-x) 



00 00 



X / dwi / dw2 / dv 



WiW2{l - x) 
Wi + W2 



v+l 



Ly{Wl,W2) ■ 



(28) 

The positivity of forward parton distributions is a conse- 
quence of the general positivity bounds on GPDs which 
will be established in the next section. On the other 
hand, we can see the positivity of the forward parton 
distribution f{x) directly from Eq. H28(l : 



00 00 



dwi I dw2 



/ W1W2 

\Wi +W2 



00 00 



dwi / dw2- 



r(j. + i) 
00 



L^{wi,W2) 



dTT" 



W1+W2 , , - 

X exp -T L^[wi,W2) 

W1W2 



X di^ [ — \ t] Li,{wi,W2). 

J \W2l — X Wil — X J 

(25) 

The step function does not vanish in the region x > |^| 
so that the above expression simplifies as follows: 



00 00 



x>\i\ 



2(1 - X)^-^ I dwi / dW2 







X dv[— + — t] L^{W1,W2). 

J \W2l ~ X Wll~ X J 

(26) 

This representation can be rewritten in the following 
form: 

00 00 

TjW (x, ^, t) = 2(1 - x)^-^ f dwi f dw2 
a:>l«l J J 





X / d'ye*'^ / diy 



r(i/ + i) 



L^{wi,W2) 



X exp 



1 1+f 1 1-e 



■W2 1 — X Wi 1 — X 



(27) 



1 



r(t. + i) 



00 oc 00 



drr'^ / dwi / dw; 



00 



X L^{wi,W2) exp ( — — ] exp ( — — ) > . (29) 

\ WiJ \ W2 J 

The positivity of the RHS follows from the inequality 

00 oc 

/ dwi / dw2Ly(wi,W2) exp ( ) exp ( ) > 

J J \ wij \ W2J 



(30) 

which is a consequence of the positivity ()19f) of the 
quadratic form L^{wi,W2)- 



VII. PROOF OF POSITIVITY 

Now we want to show that the modified double dis- 
tribution (|18|l with positive definite functions Li, (|19() 
generates GPD iJ*^^' (x,^,i) which satisfies the positiv- 
ity bounds Q). For the positivity bounds we need the 
GPD in the impact parameter representation JSJ 



rcxp 



(2^) 



■(A^6^) 



5 



(31) 



Using representation l|77|l for the GPDs i/^^' (a;, ^, i), we 
obtain 



rexp 



(2^)2 



1-e 



1 



(is exp 



47M2 



+ is(r2 - ri^ 



(37) 



Then 



1 -X 



1 



2rir2 



00 00 00 



x2(l - x)^-^ J dwi J dw2 J dj 
000 



00 00 
X J d7exp ^— J dsexp 



47M2 



X exp —7 



7" 



r(t. + i) 



X exp 



'1 



1 1+f 1 1-e 



L,, (wi, W2) ■ (32) 



W2 1 ~ a; 1 — a; 

Integrating over A^, introducing compact notation 

1 — X 1 — X 

"^1 = -r—^, r2 = 



1+r 



1-^ 



(33) 



(see Appendix^, and rescaling the integration variables 
Wk wuTk, we find 



l-x 



1 / 2rir2 \ f d'y 



00 oc 

X j dwi J dw2 exp 





2tt \ri + r2 



-i-2|^^|2_^M2(ri-r2)2 



X / dv- 



r(i. + i) 



exp 



47 



7(wi + W2) 



7-1 7-2 



rir2WiW2 



Lu {riWi,r2W2) ■ 
(34) 



Next we change the integration variable 7 — s- 7rir2 





oc 00 



X J dwi J dw2 exp 




X / az^:;T7 7T exp 



1 

47 



|6^|2-7Af2(ri -r2)^ 



(35) 



7(wi + W2) 



W1W2 



r(i. + i) 

and use the representation 

exp [-7(ri - r2)2M2] 



Lu (r-iwi, r2W2) 
(36) 



00 00 



X J dwi J dw2 exp 




is{r2 - ri) 



7(wi + W2) 



W1W2 



y. [ dvY-'"^P^,LAw^n,W2r2) . (38) 

J r(i/ + i) 

Now we turn to the positivity bound iQ written in the 
form of the integral over ri , r2 — see equation (IA8|) in 
Appendix ^ The left-hand side of this inequality is 



dri I dr2{ri + r2)^^^ p* {r2)p{ri) 







l-x 



-,N-1 



M7r3/2 ./ - y "'r(i/ + i) 





diy / d'y 



7 



1/-3/2 



exp 



(,-L|2 



47 



X ds exp — 



47Af'^ 



dri / dwi / dr2 / dw 



0000 



xLi, (wiri, ^2^2) 



P {r2) ^2 



p(ri)r, exp — zsri 



exp —isr2 

W2 



(39) 



Here we can rescale integration variables Wk — > w^/r]^. 
Then 

1 00 1 00 
dri J dwi J dr2 J dw2Li, (wiri,W2?'2) 
0000 



6 



P (Ti) r^^'^^^ exp ( —isri 



Wi 



1 



exp ~isr2 

W2 



oo oo 



dwi J dw2Ly (wi, W2) 




1 

[ , /X N+u ( ■ 77-1 

J \ Wi 

.0 

1 

/ , . N N+v ( ■ 1^2 

/ dr2P (r2 ) exp -isr2 

J \ W2 



00 00 



dwi J dw2L^ {wi,W2) (t'uiwi)(t)tiw2) > (40) 




where 



± 

(j)^{w) ^ J drp{r)r^+'' exp (-isr - . (41) 



The RHS of Eq. (|in|l is positive since (wi, W2) is posi- 
tive definite. Combining Eqs. H39|) and H4U|) . we complete 
the proof of the positivity bound (jA8|) for the GPD gen- 
erated by the double distribution lfT5|l . 



positivity and polynomiality are sensitive to the hadron 
mass M. Therefore the absence of the M dependence in 
our ansatz p8ll should mean that there must exist other 
solutions of the polynomiality and positivity constraints 
and one has to try other methods in order to find the 
other solutions. In particular, in Ref. the solutions 
of the positivity and polynomiality constraints are con- 
structed in terms of triangle perturbative diagrams. 

The parametrization of GPDs suggested here seems 
to be constructive for the model building: the positive 
definite functions L^{'Wi,W2) can be easily generated by 
using Eq. (|20|l . One should not forget about the possi- 
bility to add the D-term (|17|) which is not constrained 
by the polynomiality and positivity. 

Certainly apart from the positivity and polynomiality 
there are other theoretical and phenomenological con- 
straints on GPDs and it would be interesting whether 
representation (|18|l allows to construct viable models of 
GPDs. 
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APPENDIX A: SOLUTION OF THE POSITIVITY 
BOUNDS 

In this appendix we derive the solution @ of the posi- 
tivity bounds O . First let us define variables ri , r2 which 
can be used instead of x, ^ 



l-x 



r2 



r2 ~ n 
r2 + ri 



a: = 1 - 



2rir2 
ri + r2 



(Al) 



(A2) 



VIII. CONCLUSIONS 

In this paper we have shown that representation H18|) 
for the double distributions [understood in the sense of 
Eq. H15() ] generates GPDs (|25|l satisfying both polynomi- 
ality and positivity constraints. Our representation l|18|) 
for double distributions involves arbitrary positive defi- 
nite quadratic forms L^{wi,W2)- Functions Lu(wi,W2) 
parametrizing GPDs depend on the same amount of 
variables (wi,W2,J^) as GPDs themselves (a;,^,t). This 
means that the class of solutions of the polynomiality 
and positivity constraints found in this paper is rather 
wide. On the other hand, this set of solutions is not 
complete. Indeed, our ansatz l(TH)l does not depend on 
the mass of the hadron M whereas the positivity and 
polynomiality constraints are sensitive to M: although 
M appears neither in the polynomiality condition I^J for 
^(^) nor in the positivity bound {Tj) for F^^\ the re- 
lation © between ff^^^ and F'^) contains the hadron 
mass M. This means that the combined constraints of 



2dxd^ _ dridr2 , , 

(1 - x)^ ~ rlrl ■ ^ ' 

The region covered by the positivity bounds 

X > Id (A4) 

is mapped to the square in the ri , r2 plane 

0<r-i,r-2<l. (A5) 

Inequality |(7J takes the following form in terms of inte- 
gration variables ri,r2 (we keep variables x,^ in GPDs 
implying that they are functions of ri, r2): 

f dri [ dr2 (ri+r2\^^^ , , . , , 

— —{ P ir2)p{ri) 

J rf J \ rir2 ) 


x/'W(^x,C,^6^) >0. (A6) 



Since function p is arbitrary we can replace it 

Pin) ^ rf +3p(ri) (A7) 
which leads us to the equivalent form of inequality l)A6|l 



1 1 

dri J dr2 [ri + 




N+l 



(AlO) 



with some functions i?„. Turning back to the variables 
X, ^, we find 



Inequality (|A6|) means that function 

N+l 



ri + r2 
rir2 



(A9) 



must be a positive definite quadratic form, i.e. it has the 
following representation 



n + r2 

2^17-2 



AT+l 



(All) 

In the case of real and ^-even GPDs, functions i?„ are 
real. 

Introducing functions 



(A12) 



we obtain representation ^ for f(^) (^x,^,b^ 
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